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Abstract 

A summary is given of the work which was presented in the working group C 
"Hadronic Final States" at the DIS04 workshop (Strbske Pleso, Slovakia , April 
14-18, 2004). Progress in experimental tests and the theoretical development of 
perturbative and non-perturbative QCD is reported. 



1 Introduction 



One of the still open challenges in modern physics is to understand more completely 
the strong interaction. This is particularly important in light of the fact that many of 
the most important future discoveries in particle physics are expected to be made at 
the LHC, a machine that will collide strongly interacting particles at extremely high 
energies. 

Hadronic final states produced in high-energy particle collisions have traditionally 
been used as a testing ground for the theory of the strong interactions, quantum chro- 
modynamics (QCD). Within our working group, results on hadronic final states from 
ep, pp, NN and e"^e~ induced final states were reported by 27 experimentalists and 11 
theorists. The working group contained a vibrant atmosphere and was the forum for 
much useful discussion. 

The breadth of topics included precision tests of perturbative QCD with jet and 
leading particles. Perturbative QCD is in good shape in the regions in which it is 
applicable. However, results from HI and ZEUS on forward jets and azimuthal decor- 
relations between jets suggest that the DGLAP evolution approach may be failing and 
that calculations including small Bjorken x physics (BFKL/CCFM resummations) may 
be required. New measurements and theoretical predictions for event shapes, including 
new resummed calculations, have given improved consistency in the extraction of 
and the non-perturbative parameter ao- Amongst the main highlights of our working 
group were studies of exotic pentaquark states. In this working group, several experi- 
ments have reported their results on the search for pentaquark states with strangeness. 

2 Jets and high Pt phenomena 

2.1 Progress in experimental tests of perturbative QCD 

The production of jets and leading particles provides a clean signature for probing per- 
turbative QCD processes. At this workshop a number of measurements were presented 
from ep, pp and e~^e colliders. 

Multi-jet production in DIS provides an important testing ground for QCD. At 
HERA, data are collected over a large range of the negative four-momentum transfer 
squared, Q^, the Bjorken-x variable and the transverse energy Et of the observed jets. 
HERA dijet data can be used to gain insight into the dynamics of the parton cascade 
in low X interactions. 

ZEUS presented results on dijet and 3-jet production in DIS [1]. The data are 
described well by NLO calculations. The ratio of 3 and 2-jet production cross- sect ions 
is shown in Fig. ^ In principle, experimental and theoretical uncertainties should 
be minimised by taking this ratio. However, the theoretical uncertainty owing to the 
choice of parton distribution function (pdf) is clearly large and this requires further 
investigation. 

Insight into low x dynamics can be gained from inclusive dijet data by studying 
the behaviour of events with a small azimuthal separation between the two hardest 
jets as measured in the hadronic centre-of-mass system. Fig. |21 shows measurements 
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Figure 1: The 3/2 jet ratio in DIS as measured by ZEUS [1]. 



by HI [2] of the S distribution, i.e. the fraction of fraction of dijet events with an 
azimuthal separation of less than 120° as a function of x for different regions. NLO 
calculations for 3 and 2-jet processes are shown. Although the 3-jet calculation comes 
closest, neither calculation is able to reproduce the data at the lowest and x values. 
In this region, higher order, small x, QCD evolution effects are expected to become 
important and the validity of the DGLAP approach is in doubt.. 

A further way to experimentally probe QCD evolution is to study events containing 
high momentum jets [3] and particles [4] in the forward region i.e. next to the proton 
direction. Fig.|21shows the cross-section for forward 7r°-mesons produced at high trans- 
verse momentum as function of x. The data rise towards low x and a QCD prediction 
based on DGLAP evolution for proton structure (labelled DIR) fails to describe the 
data. Approaches based on a resolved virtual photon picture (labelled DIR-I-RES) and 
on CCFM evolution describe the data better. Analytical calculations based on the 
BFKL evolution equation also give a reasonable description of the data. 

The structure of the photon at low has been investigated by both the ZEUS [5] 
and HI [6] experiments. Fig. |3] shows the dijet cross-section as a function of for 
different regions of x^, the variable used to determine the fraction of the photon's 
momentum taking part in the hard sub-process. The data are described by the NLO 
calculations for x^ > 0.75 although for x^ < 0.75 the NLO prediction significantly 
underestimates the cross-section. This could be due to effects arising from the structure 
of the photon in this region. 

Results from the CDF [7] and DO [8,9] showed encouraging progress in analysing 
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Figure 2: The S distribution for dijet events, as measured by HI [2]. 



Run II data. The increased centre of mass energy to 1.96 GeV makes available a 
larger cross-section at the high values of jet transverse momentum. However, precision 
is severely limited due to calorimeter energy scale uncertainties for both experiments. 
The NLO predictions describe the data well although sensitivity is inhibited by the large 
energy scale uncertainty. In order to reduce the systematic effects from this source, DO 
also presented results on the angular separation between the dijets, a variable which is 
highly sensitive to the higher order emissions. This is shown in Fig. El NLO calculations 
provide an excellent description except for the extremely large (small) regions of A$*, 
where the available phase space for multi-jet emissions is limited. 

In collisions between objects with extended structure, the final state of hard parton- 
parton scattering is characterised by high transverse momentum jets but also contains 
a number of soft particles which are collectively referred to as the "underlying event" . 
This has been studied by CDF [10] in which back-to-back dijet events are analysed in 
order to isolate the rj — (p regions which are sensitive to the underlying event and to 
analyse its properties. Fig. IHl shows the charged particle density in the "MIN", "MAX" 
and average transverse region of leading jet events. On an event-by-event basis the 
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Figure 3: The forward jet cross-section, as measured by HI [4]. 



"MAX" and "MIN" regions are defined to be the regions containing the largest and 
smallest, respectively, number of charged particles. The density increases with the 
leading jet's transverse momentum in the "MAX" region and decrease in the "MIN" 
region. The complete results of this study suggest that there is a jet structure in the 
underlying event even when initial and final state radiation is strongly suppressed. The 
techniques and results of this work will be of enormous benefit at the LHC 

2.2 Theoretical developments 

Dijet data from HI and ZEUS has been a significant source of information on quantities 
such as the pdfs. In such studies it has been commonly accepted that there is a need 
to impose asymmetric cuts on the transverse energy of the two highest Ef jets to avoid 
problems of infrared enhancements. This involves cutting out a significant fraction of 
events since in fact the measured dijet total rate is largest when one allows symmetric 
cuts to be placed. Moreover, there is accurate data in the region of symmetric cuts 
and it is only the NLO QCD theoretical calculation that breaks down in this region. 
An improved theoretical estimate can be obtained by resumming large logarithms in 
the Et difference, A,^ between the two highest Et jets [11]. This allows one to use the 
theoretical calculation in the entire range where there is data. Moreover it allows 
one to test the ideas behind multijet resummations (in a problem less contaminated 
by power corrections than several event shapes) by comparing the next-to-leading 

-"^More precisely it is the difference in Et cut values. 



4 




Figure 4: The differential dijet cross-section, as measured by ZEUS. 



log (NLL) resummed theory to the experimental data. In the absence, thus far, of 
any detailed comparison of theory and data for resummations with more than two 
hard jets one would imagine that the calculation presented in [11] provides a good 
opportunity for the same. Such multijet configurations will be especially important in 
the context of hadron colliders such as the LHC As far as results are concerned it was 
shown that after resummation the slope of the total rate ex' (A) was a linear function 
of A rather than a (double) logarithmically divergent one as indicated by NLO fixed- 
order computations. As a consequence the total rate itself rises monotonically with 
decreasing A instead of turning over as indicated by NLO alone. Future extensions 
involve extending the present cone algorithm results to the inclusive kt algorithm and 
carrying out the matching to fixed-order estimates. 

Jet and dijet production in photon-photon collisions is also a valuable source of in- 
formation on the hadronic structure of the photon. While NLO calculations previously 
existed in the framework of the phase-space slicing method [12], results were presented 
here [13] for NLO computations employing the dipole subtraction method which is 
preferable, especially in the case of computing certain differential cross sections, from 
the point of view of numerical stability. The new method also allows a cross-check of 
the previous phase-space slicing results. At present the comparisons with OPAL data 
are very satisfactory for inclusive observables e.g Et spectra in single inclusive jet pro- 
duction and the direct and single resolved dominated contributions to dijet production 
see Fig. |7| However one notes a discrepancy for the double resolved dominated region. 
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Figure 5: The azimuthal decorrelation between dijets, as measured by DO. 
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Figure 6: The charged particle density for different selected regions of phase in dijet 
events, as measured by CDF [10]. 



which is still outstanding. 

Understanding the multiple collinear limit of QCD amplitudes is of great importance 
in several contexts. These include the development of subtraction terms for fixed higher 
order (NLO and beyond) computations, all order resummations, more accurate pdf 
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Figure 7: Comparison between theory [13] and OPAL data for dijet production in 77 
collisions in three different regions. 



evolution and improving the physics content of Monte Carlo event generators. For 
example the one loop correction to one to three collinear branchings is one of the 
ingredients in computing the full NNLO pdf evolution. Results were shown based 
on a new method that implements collinear factorisation directly in colour space thus 
retaining explicitly the colour matrix structure and resulting in the derivation of 1 ^ m 
splitting amplitudes. The (divergent part of) one loop splitting matrix in the limit of 
m colhnear partons reads (where si...m denotes {p\+P2 + ■■■PmY and thep^ are collinear 
parton momenta: 



5^(i)div.. _ r(i + e)r^(i-6) ifi - ( -s,j-io ' 

^"^'•••'"-^ - (4.)-r(l-2e) 2\e^,^.i^,/^ [ 



+ 



-Sl...m - 



^ i,i=l 



-7 it (7. - e^n - {la - eiT) - ^ {Po - epr 



Ki=l 



X 5p(°)(pi 



which shows the infrared and ultraviolet e pole structure. In the above the colour 
charge of the collinear parton with momentum pi is denoted by Tj, and from colour 
conservation one has J2i Ti Sp^^^ = Sp^°^ Ta {Ta is the colour charge of the parent 
parton in the collinear splitting). The flavour coefficients 7, and (3o are Iq — Iq — 
3Cf/2 and 7^ = /?o/2 = (IIC^ - 2iV/)/6. The flavour coefficients 7f ^ and p^"" are 
renormalisation scheme dependent . Sp^'^\pi, . . . ,Pm) is the corresponding tree level 
quantity. From the above equation and the tree level result one can compute one loop 
splitting functions for the 1 ^ m collinear branchings [14]. 
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3 Event shape variables 

A class of observables that are perhaps best suited for testing QCD are event shape 
variables (for a recent review see [15]). While the infrared and coUinear safety of 
these observables allows one to make perturbative predictions, event shape distribu- 
tions receive infrared and coUinear enhancements, in the form of large logarithms, that 
require resummation. Resummation tools and predictions in turn follow from a sound 
understanding of QCD dynamics in the most important kinematic regions. 

However, before one can confront resummed perturbative predictions with the ac- 
curate experimental data available for several event shapes, one usually has to account 
for power corrections. These are non-perturbative effects that, for most event shapes, 
scale as 1/Q, Q being the hard scale involved in the reaction, and make a significant 
contribution to many event shape observables. Over the past decade this problematic 
aspect of event shapes has been turned into an advantage since theoretical ideas now 
exist that permit a description of power corrections based on renormalons (for a review 
see [18]). The sizable power corrections involved in event-shape studies make this a 
valuable testing ground for ideas on non-perturbative dynamics. Hence, event shapes 
allow one to test at once both perturbative and non-perturbative aspects of QCD. 

Both the ZEUS [16] and HI [17] experiments showed new results of event shape 
distributions in DIS, measured in the current region of the Breit frame. Fig. |H1 shows 
the results of fits to event shape variable distributions, of resummed Next-to-leading 
logarithm (NLL) calculations matched to NLO calculations, with a power correction 
ansatz for the hadronisation step. The extracted values of and are shown for a 
different range of event shape variables. Good consistency for the different variables is 
obtained except for the C variable. For this variable, there is much sensitivity to the 
choice of the fit range. 

At this meeting there was a presentation on resumming Sudakov logarithms and 
power corrections for event shapes within the framework of dressed (i.e fake massive) 
gluon exponentiation (DGE) [19,20]. The basic idea of DGE is to use a massive 
gluon as the kernel for exponentiation, with the gluon mass effectively accounting for 
resummation of renormalon bubbles diagrams on each gluon: 



Here v is the event shape variable and is a Mellin variable conjugate to it while 
the suffix SDG denotes the single dressed gluon result. The above equation implies 
that the distribution for the emission of multiple independent dressed gluons can be 
obtained by 'exponentiating' the single dressed gluon result on the RHS of the above 
(in other words taking the inverse Mellin transform of e"^*-'^-'). This result preserves 
all the information present in the resummed perturbative next-to-leading log (NLL) 
results that are known for several observables. It also resums ^ (due to the dressed 
gluon), certain sub- leading logarithms whose coefficients would grow as n!, for large n, 
at n*^ perturbative order . Additionally on Borel transforming the running coupling 

^In the leading power of /3o- 
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Figure 8: Fitted values of ao and from event shape variable distributions, as mea- 
sured by ZEUS. 



as{k) in the expression for S, one generates singularities in the Borel plane, whose 
residues correspond to power corrections. 

This is contrast to the Dokshitzer- Webber (DW) [21] approach in which one ex- 
ponentiates the single massless gluon result and introduces a universal infrared finite 
extension of as to regulate the Landau singularity. In the DW approach the power 
correction then emerges as a shift of the perturbative distribution by an amount pro- 
portional to 1/Q. While the shift has been known to be an excellent practical approxi- 
mation for a wide range of values of v, it is well known that in the region of very small 
V, V ~ A/Q, the shift approximation breaks down due to increasing importance of 
higher-orders in the series of non perturbative terms A"/f"Q". The shift must then be 
replaced by a smearing of the perturbative spectrum with a shape function [22]. The 
dressed gluon exponentiation mimicks the shape function approach and, in fact, can 
be considered as a renormalon model for the full shape function. Phenomenological 
applications based on DGE can be found in [19]. For most recent results on the C 
parameter and a new class of event shapes (angularities) see [20]. 

It should be emphasised that the DGE, while attempting to maximise information 
obtained from renormalon resummed perturbation theory, is once again only a model. 
In particular it is well known that use of a massive gluon is not strictly appropriate, 
being too inclusive an approximation, for the case of event shapes. The DGE fails 
at present to address the issue of gluon decay, which is included for instance in the 
shift picture. However it is interesting to note that renormalons included via DGE 
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can reproduce the basic features of the series of power corrections resummed into non- 
perturbative shape functions. 

Regarding purely perturbative resummed results there have been significant de- 
velopments in NLL resummations for observables such as event shapes. Results were 
presented, at this conference, on automated resummations for QCD final state observ- 
ables [23] . The automation of resummation methods for generic observables in different 
processes is a big step forward since resummations done on a case-by-case basis were 
tedious and often error prone. This was in contrast to fixed order computations for 
event shapes where Monte Carlo programs exist that generate NLO results for several 
different observables without any laborious calculations required on the part of the 
user. 

The approach encoded in the program CAESAR is a hybrid analytical and numeri- 
cal method. The first step is to obtain the parametric dependence of the observable on 
final state soft and coUinear particle momenta. This must be (and for the vast majority 
of interesting observables is) of the form : 



where kt, rj and denote the transverse momentum, rapidity and azimuth wrt the hard 
(emitting) leg / and p the set of Born (hard) momenta. The parameters ai, hi and 
can be numerically obtained. 

The next step is to determine whether the observable satisfies certain conditions, 
not previously available in the literature, which guarantee the validity of standard 
approximations that underly resummation methods. The most important condition is 
that of recursive infrared and coUinear (IRC) safety which can be formulated as: 



In other words one probes the observable with two soft emissions (more precisely the 
commutator on the LHS of the above) rather than a single one (as for usual IRC 
safety) and checks that the scahng properties of the observable remain unaltered when 
one introduces the second soft emission. This condition being satisfied allows the 
exponentiation of the single gluon calculation (to NLL accuracy for global observables) 
with gluon branching merely reconstructing the running coupling in the exponent. 

The final step is to insert the parameters determined above into an analytically 
computed master formula whose form depends on the number of hard legs in the Born 
configuration (c.g n = 4 for hadron-hadron dijet event shapes). The final results are 
numerically determined by a Monte Carlo algorithm. Results were presented (see [23] 
for the transverse thrust distribution at the Tevatron collider, as an example of the 
power of the automated resummation method. 

LEP data continue to be an important source of information on event shapes as 
they are extracted in the cleanest possible environment for such studies. The LEP 
QCD working group has made important progress in quantifying and understanding 
better the theoretical uncertainties associated with event shape predictions. For details 
on recent progress made at this meeting see [24]. 




(3) 



[lim^^o, lim^/^o] -V (p, ^h, e'e^2, ■■■) = 0. 



(4) 
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4 Particle production 



The confinement transition from partons to the hadrons observed by experiments is 
poorly understood. Thus particle production phenomena cannot be explained in a 
systematic way solely from perturbative QCD. Studies of identified hadrons, inclusive 
hadronic multiplicities and particle flows allow to conduct detailed and sophisticated 
studies of many aspects of sort limit of QCD. At the same time, at this workshop we 
have seen how such observable can be used to tackle complicated problems related to 
perturbative QCD effects. 

4.1 Conventional states 

Identification of conventional baryons and mesons have received much attention in this 
working group. In particular, baryon-antibaryon production in two-photon processes 
in e+e^ collisions has stable interest at LEP community. Agreements with the diquark 
and the handbag predictions were found in reasonable shape for the L3 studies [25], 
while the expectation from the three-quark model was clearly below the data. 

Several results on the measurements of (anti)deutrons [26, 27] were presented by 
both PHENIX and HI collaborations. In heavy ion collisions, the production of 
(anti) nuclei helps to understand the production size of the system and amount of radial 
expansion when the coalescence model is used to explain the production mechanism 
of (anti)deutrons [26]. In particularly, deuteron/antideuteron spectra at mid-rapidity 
probe the late stages of relativistic heavy ion collisions. 

The measurements of (anti)deutrons are clearly a significant challenge for experi- 
mentalists, since these states can easily be produced by secondary scattering processes. 
This problem is less severe for antidcutrons which were measured for the first time in 
ep collisions by the HI collaboration [27]. The production mechanism of antidcutrons 
in ep collisions is significantly less understood compared to heavy-ion collisions. The 
result of this interesting study was in the observation that the production rate of an- 
tidcutrons in ep is by order of magnitude smaller than in heavy-ion collisions. Thus, 
in terms of the coalescence model, the size of the fireball at the thermal freeze out is 
expected to be much smaller for ep processes than for heavy-ion collisions. 

4.2 Pentaquarks 

This was the first DIS workshop at which the entire section was devoted to searches 
for pentaquarks, five-quark baryons which recently attracted tremendous attention of 
theorists and experimentalists. Let us remind that recent results from low-energy 
fixed-target experiments [28] gave evidence for the existence of a new narrow baryon 
resonance, G"*" with a mass of approximately 1530 MeV and positive strangeness. This 
state was seen in in the K^n decay channel. These results have triggered new interest 
in baryon spectroscopy since this baryon is manifestly exotic; it cannot be composed of 
three quarks, but may be explained as a bound state of five quarks, i.e. as a pentaquark, 
©"^ = uudds. A narrow baryonic resonance close to the observed mass is predicted in 
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the chiral soliton model [29]. The quantum numbers of this state also permit decays to 
K^p and K^p, which can experimentally be measured using tracking detectors [30]. 

Two DESY experiments, HERMES [31] and ZEUS [32], reported observations of a 
possible pentaquark state in the decay mode Kgp{p), see Fig. 9. The masses and widths 
agreed rather well for both measurements. However, some discrepancies with the K~^n 
decay channel [28] may exist. This ether could mean that the particle decaying to 
Kgp{p) is not the same state as that observed for the K^n channel, or, unlike the well 
reconstructed Kgp (p) final state, measurements involving detection of neutrons using 
missing mass could suffer from a significant systematic uncertainty which is not fully 
understood. 




Figure 9: Invariant-mass spectra for the Kg p (p) channel measured by HERMES [31] and 
ZEUS [32] experiments at DESY. The measurements were performed for photoproduction 
(HERMES) and DIS (ZEUS) data. 

ALEPH reported negative results on the search [33]. Presently, other LEP 
experiments have also reported their negative results (see a recent review [34]). Does 
this related to the fact that e~^e~ experiments do not have a contribution from the 
proton remnant, i.e. from the net baryon number in the initial state? At this conference 
there was a full concensus that this question requires further studies and, at present, 
no definite answer can be given. 

Note that ZEUS is the first high-energy experiment which made the observation 
of G"*" state in the central fragmentation region dominated by fragmentation of quarks 
and gluons. One should mention that, unlike fixed-target experiments, the contribution 
from the proton remnant is negligible in this region. Thus, if confirmed, this result 
needs to be understood in terms of present hadronisation models. 

It is important to point out again that several high-energy experiments have re- 
ported their negative results (see references in [34]). Of course, each negative case 
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should be considered very seriously. If the pentaquark does exist near 1520 — 1530 
MeV, possible reasons of such ZEUS success may lay in a good tracking resolution for 
K^sP (p) combinations (~ 2.5 MeV) and a relatively small combinatorial background 
for DIS at > 20 GeV^ (compared to some heavy-ion and Tevatron experiments). Of 
course, differences in the production dynamics may also be important. One possible 
way to learn about instrumental and background differences between different exper- 
iments would be to use the PDG Ac baryon (which can also decay to Kgp (p)) as a 
reference. Both HERMES and ZEUS used this state to cahbrate the mass measure- 
ments near the Kg p (p) production threshold. 

The 9+ lies at the apex of a hypothetical anti-decuplet of pentaquarks with spin 
1/2. The baryonic states Sg^g Sg^g the bottom of this antidecuplet are also 
manifestly exotic. According to the predictions of Diakonov et al. [29], the members in 
the anti-decouplct with the isospin quartet of S = —2 baryons should have a mass of 
about 2070 MeV and partial decay width into Stt of about 40 MeV. Recently, NA49 [35] 
at the CERN SPS made observation of exotic S^^"^ member of the S multiplct, with a 
mass of 1862 ± 2 MeV and width below the detector resolution of about 18 MeV. At 
this conference, several experiments, CDF, ALEPH and ZEUS reported their negative 
results on the search for '^^^^ and ^3^2. Some experiments have extremely competitive 
data with a large number of reconstructed S states used for the searches of exotic E'^^ 
and 

4.3 Multiplicities and particle flows 

The measurements of multihadronic final state allow detection of color coherence prop- 
erties. A number of interesting results were shown by DELPHI collaboration [36] using 
hadron multiplicities: a competitive value for the color-factor ratio was obtained, a 
clear signature for the destructive gluon interference was seen, and finally, the ratio 
of gluon to quark multiplicities was obtained and compared to the 3NL0 and MLLA 
calculations. Both predictions gave a reasonably good agreement to the data. 

Particle fiows were also used by DELPHI [37] to study color-reconnection effects 
between color singlets which take place in WW production. This phenomenon con- 
tributes to enhancement or depletion of soft particle multiplicity in some phase-space 
regions. The data might indicate that there is a small contribution from color re- 
connection effect, but generally, the predictions are doing rather well without extra 
contributions. A possible mass shift in the reconstructed W mass is expected to be 
rather small, by order of magnitude smaller than that expected in early days of the 
WW studies. 

Clearly, particle flows can be used to test certain details of hard QCD without 
using the jet algorithms, as long as we understand contributions from soft physics 
to such observables. There has been considerable progress in studying azimuthal 
asymmetry attributed to hard QCD by using the energy flow method as reported 
by ZEUS [38]. It was shown that this method has several advantages compared to 
conventional multiplicity-based and jet-based experimental techniques. The results in- 
dicate a rather good agreement with the NLO calculations (corrected by Monte Carlo 
models to take into account hadronisation effects). 
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4.4 Theoretical developments in particle production 

Results were also presented in this meeting, on NLO computations for inclusive and two 
photon production at hadron colliders [39]. Successful comparisons were shown of the 
NLO predictions from the program DIPHOX with data from CDF and DO Run 11 data 
for quantities such as azimuthal angle, invariant mass and the transverse momentum 
of the 7 pair in two photon production. In regions requiring resummation (such as the 
region of small pt of the pair) , as expected the comparison with NLO is less meaningful. 
However it was clear that two photon production at the Tevatron is well understood 
theoretically while a need for more data was observed in the inclusive photon case. 



5 Future measurements 

With the LHC around the corner theoretical tools have to be improved to meet the 
challenge of making new discoveries in a very complex hadronic environment. Amongst 
the most prominent tools that will be employed in such searches for new physics will 
be Monte Carlo event generators. 

Progress and developments concerning the new event generator SHERPA were pre- 
sented in this conference [40] . This is an event generator that uses the matrix element 
(ME) package AMEGIC-I— I- along with its own parton shower (PS) model and the 
CKKW procedure to merge the ME and PS. It also has an interface to the PYTHIA 
string fragmentation model for hadronization. Results from SHERPA for the jet pt 
distribution in Z+jet production at the Tevatron were shown, which indicated that the 
CKKW procedure, for such spectra, gave very similar results to those obtained from 
full NLO calculations. 

There was also a presentation on the THEPEG [41] which is a general and modular 
C++ framework for implementing event generator models. Both PYTHIA7 and HER- 
WIG++ will be built on THEPEG which will unify aspects of the event generation 
and share the administrative overheads involved. In the same talk progress was also 
reported on the merging of matrix elements and parton showers for the ARIADNE 
Monte Carlo. Results obtained in the case of the pt distribution in W+jet production 
at the Tevatron indicated that the matching procedure worked extremely well in the 
incoming quark channel and less well for the incoming gluon channel with a visible 
cut-off artifact. Work is in progress in this regard. 

Results for event shapes, jet rates, heavy quark fragmentation functions etc. were 
also displayed from the new event generator HERWIG++ [42], which was shown to 
perform as welF(and in some cases better) than the fortran HERWIG, with precision 
tuning to LEP data still to be carried out. It was emphasised that the C++ version is 
not merely a rewrite of the fortran HERWIG but there are some changes in the physics 
models. The most significant is perhaps the use of a new evolution variable q inspired 
by the heavy quark splitting function : 



ID r 2 2\ 

Pgq{z,q ,m ) = 



1 + - — 3^ 



(5) 



'Thus far only e+e processes are implemented 
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Figure 10. B fragmentation function from HERW1G++ compared to data for different values 
of the cut-off parameter 6. 



Using q as the evolution variable ensures a better coverage of the soft gluon phase 
space for emission off a heavy quark, in particular the description of the coUinear 
(dead-cone) region. Evidence for the success of the new evolution can be seen in 
the excellent description of the B fragmentation function by varying just the parton 
shower parameters (see Fig. 10). There is also an improved cluster hadronization model 
employed in HERWIG++. Extensions to include ep and especially hadron-hadron 
processes are eagerly awaited. 

6 Summary 

A selected sample of presentations given at the DIS04 workshop (working group C) 
discussed different areas of hadronic final states phenomena has been summarized. The 
detailed comparisons of the experimental data with the latest theoretical developments 
show good agreement with QCD. In spite of all these successes of QCD theory, research 
in this direction is far from complete. We hope that the increase of the collected 
luminosity at HERA and Tevatron colliders, together with future experiments at LHC 
and a linear collider, as well as further progress in theoretical QCD calculations, will 
lead to new exciting results which would open up new areas where QCD theory can be 
confronted with greater challenges. 
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